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Distribution  Results  for  Positive  Definite  Quadratic 
Forms  with  Repeated  Roots 


By 

M.E.  Bock 


SECTION  1 .  INTRODUCTION 


This  paper  examines  the  distribution  functions  of  quadratic  forms 
of  the  type 


(*) 


where  the  c^  are  distinct  positive  constants  and  the  chi-square 
variables  are  independent  with  degrees  of  freedom  Pj^  and  where 
at  least  one  of  the  p^'s  is  larger  than  one.  The  c^  are  the 
roots  and  the  p^  their  corresponding  multiplicity. 

Distributions  of  the  form  Q  arise  in  a  number  of  applica¬ 
tions.  Solomon  (1961)  noted  that  probabilities  of  hitting  targets 
frequently  reduce  to  the  distribution  of  quadratic  forms  of  the 
type  Q.  Plllai  and  Young  (1973)  show  that  the  trace  of  a  Wishart 
matrix  is  distributed  as  Q  with  all  the  p^'s  equal.  The  variable 
Q  arises  in  the  engineering  literature  described  as  a  weighted 
unbiased  Rayleigh  variate.  (See  Miller  (1975).) 

The  form  Q  also  arises  in  goodness-of-f it  tests.  Chernoff  and 
Lehmann  (1954)  showed  that  the  asymptotic  distribution  of  certain 
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chi-square  goodness-of-f it  tests  has  the  form  Q  .  The  general 
results  of  Moore  and  Spruill  (1975),  Moore  (1978)  and  LeCam,  Mahan 
and  Singh  (1983)  reveal  that  asymptotic  distributions  of  the  form  Q 
are  appropriate  for  an  even  wider  variety  of  tests.  Alvo,  Cabilio 
and  Feigen  (1982)  prove  that  the  asymptotic  distribution  of  the 
average  Kendall  tau  statistic  has  the  form  Q  . 

Section  2  of  this  paper  treats  the  general  form  Q  given  in  (*) 
with  several  p^  larger  than  one.  It  is  shown  that  its  distribution  may 
be  represented  as  a  finite  sum  of  the  distributions  of  quadratic  forms 
of  the  type  Q  with  at  most  one  p^  larger  than  one  and  with  no  more 
distinct  roots  than  the  original  form.  In  particular,  if  at  least 
two  of  the  p^  are  larger  than  one  in  the  quadratic  form  Q  ,  then 

N 

PfQ  >  c]  -  l  8  P[Q  >  c] 

.  n  n 
n=l 

where  each  quadratic  form  has  exactly  one  root  with  multiplicity 

greater  than  one.  Furthermore,  the  number  of  distinct  roots  in  each 
equals  or  is  one  more  than  the  number  of  odd  p^  in  Q  . 

A  general  method  is  given  for  reducing  the  multiplicity  p^  of 
a  root  c^  in  Q  :  For  >_  2  , 


P[Q  >  x]  -  2ClfQ(x)  +  P[Qi  >  x] 

where  fq(x)  is  the  density  of  Q  and  is  the  quadratic  form 

which  matches  Q  except  that  the  multiplicity  of  c^  is  (p^  -  2) 
rather  than  p^ 

In  Section  3  the  distribution  function  of  Q  is  completely 
specified  for  the  case  of  exactly  two  distinct  roots  c^  and  c^  , 


i.e.,  k  is  two  in  the  above  (*)  representation  of  Q.  The  distri¬ 
bution  function  is  written  in  terms  of  commonly  tabled  functions  and 
may  be  expressly  evaluated.  Section  4  relates  the  problem  of  finding 
the  distribution  of  the  trace  of  a  2  x  2  Wishart  matrix  to  the  special 
case  of  the  problem  treated  in  Section  3  where  p^  =  p^  in  the 
quadratic  form  Q  . 

In  Section  5,  results  are  given  for  the  quadratic  form 


where  0  <  c^  <  1  ,  i  ■  1,  ...  ,  r  .  This  form  arises  in  the 
asymptotic  distribution  of  certain  chi-square  goodness-of-f it  tests. 
The  distribution  function  is  completely  specified  in  terms  of  tabled 
functions  when  n  is  1  .  This  is  also  true  for  the  case  when  both 
n  and  r  equal  two  in  Qq°^  • 

The  distributional  results  obtained  in  this  work  increase  the 
importance  of  tables  for  the  distribution  of  Q  as  defined  by  (*). 
Define  P  *  p.^  in  the  definition  of  Q.  Tables  for  the  distri¬ 

bution  function  of  Q  for  P  »  2  and  3  have  been  given  by  Grad  and 
Solomon  (1955)  and  Solomon  (1960)  and  Marsaglia  (1960) .  (An  abridged 
version  appears  in  Owen  (1962).)  Johnson  and  Kotz  (1968)  give  tables 
for  P  ■  4  and  5.  Many  of  the  formulas  given  in  this  paper  reduce  the 
distribution  of  a  Q  with  a  much  larger  value  of  P  to  finite  sums 
involving  the  distribution  of  a  Q  with  a  much  smaller  value  of  P. 
Thus  these  tables  of  Solomon,  and  Grad  and  Solomon,  and  Marsaglia,  and 


Johnson  and  Kotz  become  even  more  useful.  It  seems  clear  that  the  tables 


of  Solomon  should  appear  In  a  more  easily  available  location  than  the 
present  technical  report  since  so  many  of  these  formulas  depend  on  them. 

In  the  cases  where  the  exact  finite  formulas  given  here  become 
too  involved  or  unwieldy,  it  is  clearly  appropriate  to  use  approxima¬ 
tions  to  the  distribution  function.  Solomon  and  Stephens  (1977,1978) 
have  compared  various  methods  and  give  references.  Of  particular 
interest  are  the  ones  to  Kotz,  Johnson  and  Boyd  (1967),  Imhof  (1961) 
and  Jensen  and  Solomon  (1972)  and  Davis  (1977).  The  formulas  given 
in  this  work  may  offer  insight  for  other  approximations . 


SECTION  2 


RESULTS  FOR  REPEATED  ROOTS  IN  THE  GENERAL  FORM 


Let  Cj  be  positive  constants  distinct  and  assume  that  the  Xp 
are  independent  chi-square  random  variables  with  p^  degrees  of 
freedom.  Then  define  the  quadratic  form  Q  to  be 


k  . 

I  cjXp  • 

j-1  3  j 


The  Cj  are  called  the  roots  of  the  form  Q  and  the  p^  their  multi¬ 
plicity. 

When  the  p^  are  all  even.  Box  (1954)  has  written  the  distribu¬ 
tion  of  Q  in  the  form 


k  Pj/2 


P[Q  >  c]  -  I  l  ajjSP[X2S  >  c/cj 


j-1  s-1 


'j 


k 

where  the  a.  are  positive  constants  depending  on  (c. }  and 
k  J’8 

{p  }  .  A  result  of  this  paper  shows  that  if  two  or  more  of  the 

1  i-1 

Pj  are  larger  than  one  then  the  distribution  of  Q  has  the  form 


N 


P[Q  >  c]  -  l  >  c] 

n-1 


where  each  Qr  is  a  quadratic  form  with  no  more  than  one  repeated  root 

(i.e.  at  most  one  p^  is  greater  than  one).  Furthermore  the  roots  of 

the  Q  are  the  c.'s  associated  with  the  odd  p.'s  in  Q  and 
n  j  j 

possibly  one  Cj  associated  with  an  even  p  ^ . 


Also  in  this  section  formulas  are  given  for  reducing  the  distribu- 


Vkl  kl"*"k2  y2 

Proof.  Define  W  -  I.  .  c.wx7  and  V  =  Z.A.  c...  xf  .  Then 
-  i=l  i  1  i=l  i+kg  2q^ 


W+-V  has  the  same  distribution  as  Q  and 


P[Q  >  c]  =■  P[W+V  >  c] 


P[W  >  c]  +  P[W  <  c  and  V  >  c-W] 


Now  P[W  <  c  and  V  >  c-W]  *  E^[(W<c)  •  P[V  >  c-W|w]|.  A  theorem 
of  Box  [1954]  implies  that  for  W  <  c. 


:-w|wJ. 


P[V  >  c-W|w] 


ki+k2  q± 


i-1  jjl  ai,£  *n  >  C'W^W1 


Thus 


P[W  <  c  and  V  >  c-W] 


h+k2  qi 


ii  <b  ai,iE[I(w<'=>-piV1  x“>c'w|"3 


kl+k2  qt 


1-1  Jt-X  °i.*PI\+1  X2*  ”  C"W  and  "  ‘  C’ 


Because  P[cfc  +1  >  c-W  and  W  <  c]  -  P[Q±  t  >  c]  -  P[W  >  c] 
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where 


we  have 


qi.t  •  %  +  ci+k0  ■  x2l  . 


ki+k2  qA 

>  c]  -  1  -  l  l  a  P[W  >  c] 

i=l  l-l 

ki+k2 

+  I  I  o  p[Qi  o  >  c]  . 

i-i  i- 1  ,x<  ^ 


As  noted  by  Pillai  and  Young  (1973) 


kl+k2  q± 


I  I  a  £ 

i=l  £-1  *»x 


Thus, 


kl+k2  ’l 


>  cl  ■  I  l  a  ,  ftQ,  .  >  c]  . 
1=1  l-l  X* 


Theorem  2.2.  Let  W  be  a  continuous  nonnegative  random  variable  and 
2  , 

assume  Xq  has  a  central  chi-square  (n)  distribution  independent  of 

2 

W.  Let  fg^(x)  be  the  density  of  =  W  +  CqX^  where  Cq  >  0.  For 
c  >  0,  and  n  >  2, 


P[W+c0x^  >  =1  -  <2cQ)f  (c)+P[W+c0x*_2  >  c] 


If  n  =  2,  this  is 


P[W+c0x2  >  c]  =  (2cQ)fQ  (c)+P[W  >  c] 


8 


Corollary  2.3.  For  the  quadratic  form 


Q  =  I  c  x. 

1  rt  A 


i-0 


where  the  X  are  independent  and  distinct  >  0,  i  =  0,. 


,n,  and 


Pn  >  2» 


>  c]  -  pU0)C  _2+  I  c,x!  >  c]  +  2.  -.(c) 

U  T0  i=l  1  pi 


where  fp(x)  is  the  density  of  Q  and  Xq  =  0* 

2 

Eroof  of  Theorem  2.2.  Let  ■  VH-CqX^  Let  f  (x)  be  the  density  of 


.  Then 


fQ  '<=)  -  d!  tPltffC0^  •=  c]1 
n 


We  may  write 


nw-coxj;  <  ci  - 


C/C0  n/2  -1  -u/2  fC-C0U 
u _ e 

0  r(|)2n/2  Jo 


dFTT  du 
W 


c  £*>' 
_  C  CQ 

•*0  Cc 


n/2-1 


-l/2(c-t)/c 

a 

r(|)2n/2 


0  -ft  -I 

■  J0  dF"  1 


(where  t  -  c-cpu  is  the  change  of  variable). 

Differentiating  this  last  expression  implies 


fp  (O 


■■c-t.2  (c-t)/2cQ 

—)  e 


(2C0)_1  "t I (n  —  3 )  I"  "<L--n12%  •n72.T-  [n<*wl‘ 


c-t,n/2_1  -(c-c)/2c0 

"fc  "  fc  co  6  1  fc 

+  I<n=2)  dF  +  (-1)  — -2 - 7-^n -  !  dF  dt 

Jo  _  Jo  cr(y)2n'  l_Jo  j 


(2c0)-1  (Kn>  3)P[W+c0x2_2  <  c]  +  I(n=2)P[W  <  c]  -  P[W+cQx2  <  c  J } 


Thus,  if  n  =  2 


P[W+C0X2  >  C- 


P[W  >  c]  +  (2cq) fp  (c) 


If  n  >  3, 


p[W+c0X^  >  c]  =  PfW+c0x^_2  >  c]  +  (2cQ)fp  (c) 


Corollary  2. A.  Define 


Qo  =  coxp  +  l  cii>xi 


0  i=l 


where  >  0  and  the  chi-square  variables  are  independent  and  Pq  2 


3  =  c  +  y  c^x2 

^  oxj  1£1  ci  X1 


with  density  fn  (x). 
(a)  For  odd  p^, 


P[Q0  >  c]  *  P[c^0)x2  +  J  c^l}x2  >  c]  +  (2cQ)  ^ 


(!°^) 

K  2  } 


E 

.■> 


V 


& 


\ 


J 


(b)  For  even  pn. 


£o 

2 


P[Q0  >  c] 


P[ j  c<i>xJ  >  cj  +  (2c0)  I  f„  (O 


i-1 


j-1  2j 


Proof .  A  recursive  application  of  Corollary  2.3  gives  the  result.  | 
Substitute  the  result  of  Corollary  2.4  in  Theorem  2.1  to  obtain 


the  next  result,  recalling  that  ai  £  * 


Corollary  2.5.  Define  Q,Qq,cx^  ^  and  c^  as  in  Theorem  2.1.  Then  if 


fi  ^(y)  is  the  density  for  the  quadratic  form  Qo  +  ci+k  *2£’ 


we  have 


kl+k2  qi  l 

P[Q  >  c]  -  P[Q0  >  c]  +  l  l  B  f  l  f  .(c) 
U  i-1  £-1  j-1  l»J 


where  6i>£  -  2ci+kn  ai,JT 


Theorem  2.6.  Let  {c.}n  be  distinct  and  positive  and  let  the  chi-square 
-  i  1=0 

variables  be  independent.  For  c  >  0  and  Cq  >  c^,  i  *  l,...,n. 


n 


P[  l  c  x:  +  C  X,  >  C]  -  P[  l  CX  >  c]  + 
i-1  1  Ti  U  i=l  1  Ti 

+  PtI  cV  <  c] 
i-1  Ti 


o  ;  r^2 

i-1  lCH 


-1 


where  c^  *  (c^-Cq^) 


11 


SECTION  3 


QUADRATIC  FORMS  WITH  TWO  DISTINCT  ROOTS 

In  this  section  we  consider  quadratic  forms 

Q  “  Cl^  +  c2^n 

where  the  c^  >  0  are  distinct  and  the  chi-square  variables  are 
independent. 

The  density  of  Q  is  given  in  Theorem  3.1  in  terms  of  a  confluent 
hypergeometrlc  function.  (It  is  given  in  Miller  (1975,  page  61)  in  essen¬ 
tially  the  same  form.)  Simple  representations  for  the  confluent  hyper¬ 
geometric  function  whose  first  two  arguments  are  integers  divided  by 
two  are  given  in  Bock,  Judge  and  Yancy  (1984)  and  may  be  used  to  evaluate 
the  density. 

Finite  representations  for  P[Q  >  c]  are  given  in  terms  of  commonly 
tabled  special  functions.  In  the  case  that  both  m  and  n  are  even, 
the  result  of  Box  (1954)  given  in  section  2  shows  that  the  special 
functions  are  probabilities  of  chi-square  variables  with  even  degrees 
of  freedom.  (There  are  simple  formulas  for  these  given  in  terms  of  the 
exponential  function  and  a  polynomial.)  Theorem  3.3  provides  a  finite 
representation  when  m  is  odd  and  n  is  even.  If  c ^  >  c^»  the  special 
functions  are  probabilities  of  chi-square  variables  with  odd  degrees  of 
freedom.  If  c^  >  the  special  function  is  Dawson’s  integral  tabled 
in  Abramowitz  and  Stegun  (1964).  For  odd  m  and  n.  Theorem  3.5  implies 


that  P[Q  >  c]  may  be  written  as  Pfc^x^+c^Xi  >  c]  (tabled  by 
Solomon  (I960))  plus  a  finite  number  of  modified  Bessel  functions 
(tabled  by  Abramowitz  and  Stegun  (1964)).  Furthermore,  Ruben  (1960) 
notes  the  (unpublished)  result  of  Dr.  David  Kleinecke  that 


P[c'1)xM2>XJ<cl 


Pt*2,A<B1  -  PtX2,B<Al 


where 


C,  -1/2  r~1/2\2  R  -  C,  -1/2.  -1/2  2 
4(cl  2  ’  *  B  4(cl  ^2  ' 


and  X2  \  is  a  noncentral  chi-square  random  variable  with  two  degrees 
of  freedom.  Extensive  tables  of  noncentral  chi-square  probabilities  by 
Haynam,  Govlndarajulu,  Leone  and  Siefert  (1983)  and  these  may  be  used  to 
supplement  the  tables  of  Solomon  (1960)  and  Marsaglia  (1960) . 

Theorem  3.1.  Let  m  and  n  be  positive  integers  and  let  c^»c2  and 
c  be  positive.  Then  the  density  of 

Q  "  Cl4  +  C2^ 


is 


fn<y> 


(m+n)/2-l  -y/2c. 


r(^)(2c1)m/2(2c2)n/2 


lFl(f,  ^(c^-c"1)  \) 


2  2 

for  y  ^  0  where  and  ^  are  independent  chi-square  random 


variables. 


Proof.  Let  and  W2  be  independent  random  variables  such  that 

W^/c^  ^as  a  c^“scluare  (m)  distribution  and  W2^c2  ^as  a  cM“S{luare 
(n)  distribution.  Then  the  density  of  is 


h1(x) 


m/2-1  -x/2c 
x _ e _ 

r(§)  (2c1>m/2 


1 


for  x  ^  0. 

The  density  of  is 


h2(x) 


n/2-1  -x/2c2 
x  e 


r(|)(2c2)n/2 


for  x  _>  0. 

Then  the  density  of  Q  *  W.j+W2 


fQ(y) 


7 

h. (y-x)h_(x)dx 

h 


-y/2c. 


I): 


y  ,  .m/2-1  n/2-1  -‘/^(n^-n-1) 

(y-x)  x  e 

0 


(2Cl)m/2(2c2)n/2r(|)r(|) 


The  integral  in  parentheses  can  be  written  as 


$n+n)/2-l 


r(f)r(|>y 

r(S^) 


1  ll2’  2  ’ 


according  to  Lemma  1  of  the  Appendix.  Thus  for  y  ^  0, 
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fQ(y) 


-y/2c.  (m+n)/2-l  , 

e  y  1F1  (c^  -C2  )  2") 

r(2±H)(2c1)m/2(2c2)n/2 


The  following  is  a  direct  result  of  Corollary  2.3  and  Theorem  3.1. 

Corollary  3.2.  Let  Cq,c^  and  c  be  positive  and  assume  that  chi-square 
variables  are  independent  in  the  following  expressions.  Then  if  m  >  2, 


?lcoxl+ci>i  >  c]  "  p[coV2+cA  >  c]  + 


(nH-n)/2-l  -c/2Cq 

C0(2} _ 1 _  p  (n  HtS..  c,  -1  -1... 

r,BHu  m/2  n/2  1  P2’  2  ’  2^  0  1  ’’ 

r(“)co  c0 


For  m  ■  2, 


Ptc0*2+ClXn  >  C1  *  >  fH  + 


c  0/2 

V2c/  -c/2cft  .  . 

1  „  0  /H  n  ,  ,  c  f  — 1  — 1* . 

r<|+i)  112'2  ’  2  0  'ci  ' 


Thereom  3.3.  Let  m  be  an  odd  integer  and  let  n  be  an  even  integer. 
Then  for  positive  Cq  and  c^  and  c,  and  independent  chi-square 
variables, 


P[cOxm+clxii  >  *  Pf>U  =•  +  <*> 


c  -1  -1 

If  c^  >  Cq,  then  x  m  2  ”co  )  <  0  and 
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(*) 


r(f) 


rtt  +  ?  -.  2 


m/2  n/2-1  ^  - ±_  pfv^  >  — ^— ] 

JQ  P‘Xm+22,  <  _2xJ  L’x2  Xn_2Jl  2cl 


(The  above  also  holds  if  m  is  even.)  If  <  Cg,  then 
x  -  |  (e^-c'1)  >  0  and 


(*)  *  e 


T^.2t  >  ■%  1  •  at  •  -  lHT'2 

il-o  ^  C1  l  /frM  t-0  r(t+y) 


where  V(y)  is  Dawson's  integral  and 


(-l)(m“1)/2r(!+ Jl)  m/2 

- - - - i -  (^_)^(^) 

l  r(®)£.  2cix  2cqx 


Proof.  Let  m  be  an  odd  integer  and  n  an  even  integer. 


P[C0^+C1X!  >  C] 


-  > f;i  +  ^ <  %  a”d  4  >  !Z^a' 


The  last  probability  in  the  sum  is 


(*) 


c/c_  m/2-1  -u/2 


_  fc/c0 

Jo  r(|)2 


m/2 


°°  n/2-1  -x/2, 

x _ e  dx 

(c-c0u)/Cl  r(^)2n/2 


}  du 


The  inner  integral  above  in  brackets  is 


-> 


'■'.'*  '  .'  7 


Using  this  and  making  the  change  of  variable  w  *  cqu/c  in  (*) 


implies  that 


-c/2c. 


(*) 


s - (-£-)  y 

r<?)  2V  jio 


(— ) 

m/2  n/2-1  ^2c/  rl 


j!  J 


m/2-1,,  v j  -"c<c01”cl1)/2 
w  (1-w)  e  dw 


*c/2cl  m/2  n/2-1  v2c 

e -  /_£_)  y 

ref)  2co  j-o 


(^)j 


7—  [B(f,j+l)1F1(f,J+l+f;  f(c‘1-c01))] 


The  evaluation  of  the  integral  follows  from  equation  3.393,  #1,  p.  318, 
of  Gradshteyn  and  Ryzhik.  Thus 


(*) 


-c/2c,  m/2 

e  1(— ) 

zCq  n/2-1 


c  J  r(f> 


r(f) 


l 

j-o  2ci  r(f+j+i) 


T-i  Jil  I  ^  C  /  "1  “1\  \ 

1F1(2,:5+1  + 2’2(cl  ~C0 


c  -1  -1 

If  Cq  >  c^,  then  x  *  ^(c^  -c^  )  <  0  and  Case  2'  of  Theorem  2  (in 
Remark)  of  Bock,  Judge  and  Yancey  (1984)  Implies  that 

m  (|,j+l+|;x)  (_1)(m-l)/2  j-Km+l)/2  T(k -•|)x_(k'1)P[X2k_1  < -2x] 

2>1  1  T(|+j  +  l)  (-x)1/2  k- (m+l ) / 2  (k  -  -^±±l):(j _  k). 


Interchanging  the  summations  in  (*)  we  may  write  for  l  -  k- 


(m+1) 


<*)  - 


'C/2ci  .  _/2  n/2-1  T(Z  +  f) 

2 _  /  -c  -.m/2  y  _ 2 

T(^)  2c0X  i-0  ZlxZ 
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^\n+2£ 


<  -2x] 


n/2-1  (2c  } 

l  1 


J-* 


(J-*)! 


then 


j  z!_ 

n/2-1  ^2c  ;  e  1  2c 

Not.  that  ^  -jpjyy  la  Hv2,  >  I£  '0>C1> 

c  -1  -1 

x  -  -j(c^  ~cq  )  >  0  and  by  Case  2  of  Theorem  2  of  Bock,  Judge  and  Yancey 
(1984) ,  implies 


T(|+ j+1) 


^(f.J+l+fix)  -  e*(-l)<^1>/2 


t  ,-1/I‘M!  tlltf.,)  ^1/2,  (»-3)/2+J-s  (.x)t  , 

*  L.  s!  H-s^!  1  r—  L  -x  _* 


(j_s):  /!^  t-o  r(t+4) 


(  n(m-l)/2  \  (m-l)/2  r(&  +  2}  2P(x1/2)  &+(m;3)/2  (_x)t 

('1)  Jo  /s  ‘  Jo  r(t+f> 


where  V(y)  is  Dawson's  integral,  tabled  in  Abramowitz  and  Stegun. 


.'c/2c°  <-£> 
Sfr  u 

o 


(nH-l)/2 


«/2“l  c  ^  i 

(-D  I  (^")  I 

j-0  *C1  4-0 


(m-l)/2 

Wf) 


C— N  4-n  2c.  o-n  4!(j-4)I 


1/2,  4+(m-3)/2 

t  )  _  y  .  (-x) 

*  t-o  r(t+|) 


-c/2c 


0,  c  ,m/2 


2C°  ,ZP(*1/2>  {n£2-1  _c_  j  1  HUf) 

(-l)(m_1)/2r(y)  /nx  j-0  2cl  4-0  4!(j-4)!x£+(m_1)/2 


n/2-1  .  1 

-<  I  (^-)j  I 


r(4+|) 


4+(m-3)/2  ,  tt 


j-o  2ci  juo  4:(j-4):xU(ra_1)/2  t-o  r(t+4> 
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By  interchanging  the  orders  of  summation  the  subtracted  bracketed  expres¬ 


sion  with  three  summations  may  be  written  as 


/_c_\  ^T(o 

n/2-1  4^'  ^  2'  H+( m-3)/2 

JUO  A!xi+(m-1)/2  t-0 


j-£ 

n/2-1  v2c1; 

ih 


r(t+f) 


sc/2cl  n/2-1  r(£+|)  c 
xu/2  V.  (2cr 


+  !>  ,  c  £+<m;3>/2  (-!)< 
ll  ''2c,x)  i 


t+1/2 


r(t+|) 


*  P[ Cll  >  * 


A  similar  Interchange  in  the  bracketed  expression  multiplying  V(x1/Z) 


implies 


x/2(  .(m-l)/2  c,  m/2  «,  1/2.  f n/2-1 

<*>  “  6  W -  (-i)  J  J 

r(f)  c2  ^74  1  1-0 


ra+f)  r 

ll  '■2c1x; 


H+m/2 


P[C2£>^f  + 
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n/2-1  r(£  +  S)  £+m/2 

Y  _ ±—  (-£ _ ) 

£-0  ''2c1x; 


P^-2£  >  ^ 


£+(m^3)/2  (j>1)txt+l/2 

t-o  r<t+|) 


n/2-1 


<*>  -  •*  ’  l  nCii  >  t]  ’  “« 


1p(,1/2)  .  *+<"j3>/2 
.  /F M  t-o  r(t+-|) 


where 


r(f+£)(-l)(m"1)/2  A  m/2 

— - 5 - Cjf-j)  Czf-j)  •  (I 

r(|)  :  2cix  2co 


Remark:  (Following  Theorem  3.3)  Rummer's  transform  (^F^(a,b;c) 

e  i^iO>-atb;-c))  implies  that  (*)  may  be  written  as 


-c/2c. 


(*)  “ 


r(f+i) 


,  C  v.1 

«  c  m/2  j  n/2-1  ^c^  r(j+l)r(l  +  |) 

:“+i)  <2co)  |  j-o  j:  r(i)r(j+i+|) 


Note  that  the  summation  in  the  bracket  represents  the  first  ^  terms 

in  c 

of  an  infinite  series  expansion  for  jF^d.l  +  y;  )  given  by  Slater 
(1960)  equation  (2.3.2),  p.  22.  Thus  as  n  -*■  00 ,  the  sum  becomes 

lF1d»l  +  2;  Icjj*’ 
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0 


(a)  For  cx  >  cQ, 

Ptc0^+clX2  >  ^  “  pt^  >  ^  + 


e 


-) 

0 


m/2 


<  cCc^-c-1)]  . 


(b)  For  c1  >  Cq t 


Pl<:0V+‘=lX4  >  '1  '  P^  >  + 


-c/2c 


c.  m/2  a  -t-i 

{P^<e(c-1-c-1)](l  +  ^-)  + 


-1  -lxl  m 


P'C2 «  c(co  -*1  51  f 


}  . 


(c)  For  Cj^  <  cQ  and  x  -  j  (c^-c*1). 


22  2c  _c/2co  -c 

^0W2  >  c]  -  P[^  >  ^]  +  e  °(^) 


Jp(x1/2)  _  (m^3)/2  (-pVl 

1  /ro74  t-o  r(t+|)J 


Corollary  3.4.  For  Independent  chi-square  variables  and  positive 
and  c,  and  x  -  y  (c^-Cq1) 


(a)  if  c,  >  c_  • 


then  x  <  0  and 


Y,v| 


Pfc0Xl+ClX2k  >  c]  *  P[X1  >  ^  + 


-C  1/2  ky 1  P[xi+2&  <  ~2x]  2 

(2vx  alo  a  1  -1  Xk“2£  c  1 

u  560  i!xVa  +  |))  1 


(b)  if  <  Cq,  then  for 


.  r(&  +  2)  ,  c  .A,  c  .1/2 

l  "  V.  ^CjX^  ^CqTTX 


P ^ CQXl+ClX2k  >  c^PCxJ  >  fH  + 


**  ^  P^X2(k-£)  >  “1  *  «A 

SL-0  C1  *  /f/4 


Jl-1  .  t 
£  (-1)  x 

t=o  r(t+|) 


Examples  (for  Corollary  3.4).  Assume  <  c^.  Then 


P[c0Xl+ClX2  >  c]  -  P [xj  >  “]  + 


-c/2c-  4c,  1/2 


+  •  p«f  <‘lS  »  > 

O  1 


p[coxi+cix4  >  c]  *  p[xi  >  + 


-c/2cn  1/2  , 

+  6  (2^}  C2P(x1/2>  (1  +  ~(l+x  X)) 


2c1  (x 


where  x  *  j  (c^-c”1). 


If  both  m  and  n  are  odd,  then  Theorem  2.1  implies  that 


(m-1) / 2  ,  2 

pEQ>c]«  altl  p[CiXi+2£+c2Xi  >  C1  + 


(n-l)/2 


2.  2 


+  2=1  a2'1  P[ClXl+C2Xl+2£ 


where 


c,  (m-1) / 2  c.  (n-1) / 2 

a  .  =  ( — — )  ,  a  =  ( - — ) 

1,  C1  c2  2  n~^~  c2-ci 


(n-1)  ,  c2  /  „  ,  m-3 

^  ml  '  a  ml*  OO  ^  ^  "  !»•••»  > 

1,^-i  1,  2=1  2£  C2“C1  2 


(m-1)  ,  C1  /  0  ,  n-3 

2  ml  t  »  nzi  *  2i  *  c  -c  ,£  2  ' 

*9  2  ^  ^  >  2  i.  ^ 


The  next  theorem  evaluates  the  probabilities  in  the  summations. 

Theorem  3.5.  Let  m  be  an  odd  integer  larger  than  one.  For  independent 
chi-square  variables  and  positive  c^,  and  c, 

p[co^+cixi  >  c]  =  p[co°)xi+ci1)xi  >  C]  + 

-(d0+d1)(m-l)/2  d-ij!  i  k 

6  jk  r(j+I)  iio  (2'k)!  do 

where  d.  =  c/ (4c . ) ,  i  =  0,1,  and  I  (t)  is  the  modified  Bessel  function 
i  l  n 

Remark:  Note  that  1  (t)  =  I  (t)  and  1  (-t)  =  (— 1) nI  (t).  Using  the 

-  -n  n  n  n 


recursion  formula 


m  LiUi-U 


|W| 

\ 


I  ! (t)  =  I., (t)  +  f2-  I  (t) 
n-1  n+I  t  n 


all  the  Bessel  functions  may  be  written  in  terms  of  I^Cd^-d^)  and 

VW 


Proof.  Corollary  2.4  implies  that 


(m-3)/2 

>  c]  -  p[c0°)xj+c{1)xj  >  CJ  +  (2c0)  lQ  fQ3+2.(c) 


.(0)  2  n>  2 


where  f  (c)  is  the  density  of 


Qk  C0Xk  +  C1X1 


Theorem  3.1  implies  that 


3-+2j+l,_  x  _c/2c 

2  1  c/  c0  _  .1  3+21+1.  ,  -1  -he. 

c  e  2  *  (co  ^ 2^ 

%+21  ./2 
J  r(^P)(2c0)  2  (2c,)1' 


cj+le  c/2c0  ^<1  J+2i  (c^-c-1)!) 


(j+l)!(2c0)j+3/2(2c1)1/2 


Kummer's  transform  for  confluent  hypergeometric  functions  implies 


y 

b 

$ 


fQ  (O 
W3+2j 


e  cJ  1F1  ( ( j+1)  +-j,  (j+l)+l,  (cx  -c( 
(j+1) !  (2cQ)  (j+D+1/2(2Ci)l/2 


(m-3)/2 

Thus  (2c_)  l  f_  (c)  is 
0  j*0  Q3+2j 


-c/2c.  c  1/2  (m-l)/2  (2cn}  . 

2  )  I  "XT—  iFi(k+-|,k+l;(c 

1  k«l  k.  i  x  z 


Lemma  2  of  Che  Appendix  implies 


(m-3)  /2  ire 

(2c  )  l  f  (c)  =  (--£) 
0  j=Q  Q3+2j  C1 


irc0  1/2  -c/4(c“1+c11) 


(m-l)/2  (4cn)  k! 


l)/2  '*c„'  -  k  1  la^c-1- 

k=l  r(k+i)  j=0  (k_3^!  1=0  £!(j-«.)! 


Adding  this  to  P[Cq0\2+c^x2  >  c]  gives 


P[c0^+ClXl  >  C]  for  di  =  1  =  °‘1 


Examples  (for  Theorem  3.5). 

Set  dA  -  c/(4ci),  i  -  0,1. 

(a)  p[c0X3+ClXl  >  c]  -  P[c(J1)x^+c^2)xJ  >  c]  + 


_ _  -(dn+d  ) 

+  ’^77 e  (Io(drdo)+Ii(drdo)) 


pIcOX5+cixl  >  c]  =  ptco  xi+ci2)xl  >  c]  + 


, _  -(W 

+  ■HdI e  + 


+  Ii(drdo)[1+l  do!  +  I2(drdo)[32') 


>  c]  +  4  e  W°  dl) 


U  (d  -d  Hl+f  d0]  +  I1<d  -d  )[i+-j2  ( 


SECTION  4 


THE  DISTRIBUTION  OF  THE  TRACE  OF  A 
2-DIMENSIONAL  WISHART  MATRIX 

The  distribution  of  the  trace  of  a  Wishart  matrix  arises  naturally 
in  problems  that  involve  a  normal  sample  covariance  matrix.  Let 
i  ■  1,...,N  be  a  random  sample  of  p-dimensional  vectors  from  the 
normal  population  N(y,E).  For  X  m  E^^  X^N,  it  is  well  known 
that  the  sample  covariance  matrix 

N  _  t  _! 

S  -  l  (X .  -X)  (X  -X)  C  (N-l) 
j=l  2  2 

has  a  Wishart  distribution  (N-l, (N-l)  ^E).  (See  Anderson  (1958).) 
Pillai  and  Young  (1973)  have  shown  that  tr(S)  has  the  distribution 


where  the  are  the  characteristic  roots  of  the  positive  definite 

matrix  (N-l)  and  the  ^  are  independent  chi-squared  variables 

each  with  N-l  degrees  of  freedom. 

In  this  section,  the  distribution  of  tr(S)  is  analyzed  for  p  =  2. 

Theorem  4,1.  Let  the  chi-square  variables  be  independent  and  c^,  C2 
and  c  positive.  For  an  odd  integer  m  larger  than  one, 


<>*««>**  > 


2.‘<‘‘1+d2)(M1d,)(«-«/2 

+  1  d  d  («-3)/2  (dl'd2) 

tr(f))‘  rwi^)  <2> 

[d  +d,] 

+  [d^I  1  (m~l)  (dl_d2)}’ 


where  di  -  c/Uc^,  i  -  1,2. 


Proof .  Two  applications  of  Theorem  3.2  (first  to  and  then  to  c2) 

yield 


Ptcf >  c|  -  Plcl<1)^-2+<=i2>>4.2  >  =1 


m-1  -c/2c. 

/Cv  X 

C1  'T'  e  m  11^ 

+ -  lFl(2,m;(cl  _c2  )f) 

r(m)(Clc2)m/2 


.  _  o  -c/2c 
/Cvin—  L  1 

®  _  _i  _i 

»  *  t?  /5L  i  i  c  /  “X  ”  x v , 

„  v (m-2)/2  1Fl(2‘1,m_1;2(c2  _C1 

1  (nt;  (c^c^  c2 


Setting  d^  »  c/ (4c^) ,  i  =  1,2,  we  may  write 


<*>  •  Pt«=i1)X2+=22)>£  >  “=1  -  PI^1)W+d22)>W-2  >  d> 


(2d1)m/2"1(2d2)m/2  -2d 

- -  e  JJ^.m^Cd,^)) 


(2d1)m/2_1(2d2)m/2“1 

foS^T) 


-2d 


2  1F1(|-l,m-l;2(d2-d1)) 


Rummer's  transform  Implies 


-2d„  -2d 

*  Vl(f'1’m_1’2(d2'dl))  "  e  1F1(|,m-l;2(d1-d2)) 


-2d,  (d-d-) 

“  e  t1F1(f-l,m-2;2(d1-d2))  +  ^(f.m^dj-d^)} 


by  Lemma  3  of  the  Appendix.  Substitute  this  in  (*)  and  write 


(*) 


-2dl 

e  .m/2-1..-,  .m/2-1 

rlmT  (2di)  (2d2> 


{(2d2)  1F1(|,m;2(d1-d2)) +  (m-D^F^S-  l,m-2;2(d1-d2>) 


<drd2> 


_2dl,/.  .  .m/2-1 
e  (4djd2) 


rfej -  {Wl+d2>lFl<!-"!2(drd2)) 


Equation  12.6.3,  p.  309  of  Abramowitz  and  Stegun  implies  that 


^(^^(dj-dj)) 


wSti,  drd2drd2  -(-1)/2  T  , ,  ,  , 

r(  2  )e  (  2  )  I.m-1. (dl~d2) 

(  2  ’ 


and 


1F1(|-l,m-2;2(d1-d2))  -  rPj^e 
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m-l^drd2  dl_d2  “(tn“3)/2I  (<*,-<*2 

(  2  }  (Si3)  1  2 


LlaJ-a  r-j 


Substituting  these  in  (*)  we  have 


(*) 


■(dl+d2)rrm+lw/  .m/2-1  ..  . 

e  r(— =— )(4d  d  )  (d.+d-) 

- 2  1  2 -  2{tM  I 


d.-d,  m-3/2 

r(m)(-~ 


(d,-d0)  +  I  _  ,  (d,-d0)} 


(drd2)  (^>  1  2'  (^r1  2' 


Corollary  4.2.  Let  m  be  an  odd  positive  integer  larger  than  one.  If 
c^,  c2  and  c  are  positive  and  the  chi-square  variables  are  independent 


then 


He^xJ+c'2^  >  01  ’  P[ci1)Xl+cl2)Xl  >  '1 


-(d1+d2) 


m-1 


2(d1d2)1/2  <2J> 


2  ,,  8d.d  j 

l  T2jtr  (d^d!>  Udrd2)i  (dl-d2)+(d1+d2)i  (dl-d2)} 


‘1  2 


where 


d±  -  c/(4Ci),  i  -  1,2. 


Examples  (for  Corollary  4.2). 

(a)  P[cJ1)X2+cf  }x2  >  c]  -  P[c^1)x2+c^2)X2  >  c] 


+e 


-(d1+d2) 


<Mid2>l/2(Io(drd2>  +  VW1 


(d1+d2) 


(b)  P[c^\2+c^2)X2  >  c]  -  P[cJ1)xJ+c^2)xJ  >  c] 


-(d  +d2)  ,  I.Cd.-d,)  o 

+  «  <4dld2>  (l0(drd2>  +  -(■d1--d2)-’dl4d24I  dld2> 
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•>-' 


A'. 


SECTION  5 


APPLICATIONS  TO  CHI-SQUARE  TESTS 


When  sample  maximum  likelihood  estimates  are  used  in  chi-square 
tests  for  goodness  of  fit,  Chernoff  and  Lehmann  (1954)  have  shown  that 
the  asymptotic  distribution  of  the  test  has  the  form 


* 

Q 


(i)„2 


i-1 


+  XP 


where  the  c^  satisfy  0  <  c^  <  1  and  are  distinct  and  the  chi-square 
variables  are  independent.  For  p  =  2,  Corollary  5.2  gives  the  exact 
percentage  points  of  Q  in  terms  of  the  distribution  of  c'^X^ 

which  may  be  taken  from  the  tables  of  Solomon  when  r  =  2  or  3  or 
the  tables  of  Johnson  and  Kotz  when  r  =>  4  or  5.  See  also  the  formula 
of  Kleinecke  given  in  the  first  part  of  Section  3.  For  p  =  4,  and 
r  ■  2,  Corollary  5.4  gives  the  distribution  of  Q  in  terms  of  Solomon's 
tables  and  modified  Bessel  functions.  When  p  ■  4  and  r  =  3,  Corollary 
5.3  gives  the  distribution  of  Q  in  terms  of  Solomon's  tables  and 
Johnson  and  Kotz's  tables.  Moore  and  Spruill  (1975)  give  similar  limiting 
distributions  for  general  chi-squared  tests  of  fit  as  well  as  a  unified 
approach  to  the  problem.  An  example  of  Moore  (1978)  for  a  chi-square 
test  of  fit  with  random  cell  boundaries  and  M  cells  gives  a  limiting 
distribution  of  the  form 

Q  “  X(M-2)  +  Xxl  * 

It  is  clear  that  if  M  is  even,  Corollary  3.4  allows  one  to  compute 


exact  distribution  points  of  Q  from  standard  chi-square  tables.  If 
M  is  even  Corollary  3.5  implies  that  the  use  of  Solomon's  tables  and 
those  for  the  modified  Bessel  functions  allow  the  exact  computation  of 
percentage  points.  (If  M  is  large,  it  is  better  to  use  an  approxi¬ 
mation.  ) 

The  asymptotic  distribution  of  (l+nTn)  where  ~n  is  the 
average  Kendall  tau  statistic  has  been  shown  by  Alvo,  Cabilio  and 
Feigen  (1982)  to  be  of  the  form 


ClXr-l  +  c2Xm 


where  m  ■  (r_^)  end  —  *  r+1  and  the  chi-square  variables  are  inde- 

C2 


pendent . 

See  the  results  of  Le  Cam,  Mahan  and  Singh  (1983)  for  other  appli 
cations. 


Theorem  5.1.  Assume  Cq  >  c^  >  0,  i  =  l,...,n,  and  the  chi-square 


variables  are  Independent.  Then 


P[  j,  c<»xKxL  »  '  P‘  j,  >  '1 


i-1 


i=l 


i=l 


j*l  k=0  0 


(k  }n 

1  £JJ.-1 


n 

I  ko  =k 
£=1 


0  <  k?  <  k 


n  2c*  ki  T(k  +h  n 

(n  (--f)  — >p[  l  c  x  v- 

i-1  c  k . !  r  (4)  i=i  1  f 


i,  -  •  "i  ■  2'  r  .*..2  <  c]\ 


and  K  =  e 


where  ci  -  (cjL  -cQ  ) 


n 


i=l  0 


Proof. 


PI JiCi1)xl+C°X2»  *  Cl  '  PIJ1^1>^  '  ^ 


n 


(*)  =  P[  l  c^X2  <  c  and  c~X 
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By  Lemma  4  of  the  Appendix,  for  W'  =  c  X^  and  c 
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We  can  note  that  the  desired  result  is  true  since 
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for  any  real-valued  function  h,  i.e. 
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This  following  corollary  follows  from  Theorem  2.6  or  Theorem  5.1. 


Corollary  5.2.  Let  c^  >  c^  >  0,  i  =  l,...,n  and  the  chi-square 


variables  are  independent.  Then 
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j  c  -  (c  -cn  )  and  K  =  e  U(  II  (1 - -)) 

110  i-1  CQ 

The  following  corollary  follows  directly  from  Theorem  5.1  for  m  = 
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Corollary  5.3.  For  Cq  >  >  0,  i  =  l,...,n  and  independent  chi-square 

variables,  for 
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In  the  following  corollary,  Solomon's  tables  and  the  modified  Bessel 
function  tables  allow  one  to  evaluate  the  distribution  of  the  quadratic 
form. 


Corollary  5.4.  For  c^  >  c^  >  c2  >  0 
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Lemma  1. 


fU  v~1/  \P-1  Bxj 

x  (u-x)  e  dx 

J0 

„  U+v-1  r (y)T(v) 


F, (v,y+v;8u)  for  y,v  positive. 


r (y+v)  lrl 


(See  Gradshteyn  and  Ryzhik  (1965)  equation  3.383,  //l,  p.  318.) 


Lemma  2.  For  real-valued  c  and  positive  integer  n. 
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Proof.  Equation  (13.4.9),  p.  507,  of  Abramowitz  and  Stegun  (1964) 
Implies  that 
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The  product  theorem  for  differentiation  implies  that 
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Equation  9.6.29  of  Abramowitz  and  Stegun,  p.  376,  is 
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l.  Let  m  be  an  odd  positive  integer  larger  than  one.  Then 
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